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Abstract 

We study a generic class of models for ultra-high energy cosmic ray (UHECR) phenomenology, in which the 
sources accelerate protons and nuclei with a power-law spectrum having the same index, but with different values for 
the maximum proton energies, distributed according to a power-law. We show that, for energies sufficiently lower than 
the maximum proton energy, such models are equivalent to single-type source models, with a larger effective power 
law index and a heavier composition at the source. We calculate the resulting enhancement of the abundance of nuclei, 
and find typical values of a factor 2-10 for Fe nuclei. At the highest energies, the heavy nuclei enhancement ratios 
become larger, and the granularity of the sources must also be taken into account. We conclude that the effect of a 
distribution of maximum energies among sources must be considered in order to understand both the energy spectrum 
and the composition of UHECRs, as measured on Earth. 
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1. Introduction 

The average nuclear composition of ultra high energy 
cosmic rays (UHECR) at their sources is one of the key 
ingredients in their phenomenology. It has been shown 
that a mixed composition, that is to say the presence 
of nuclei in addition to protons, implies a harder source 
spectrum, typically in E~ x with x 2.2-2.3 as opposed 
to x a 2.6 - 2.7 in the case of a pure proton composition 
(depending on the cosmological evolution of the source 
power) H). In addition, a mixed composition implies 
that the energy at which the extra-galactic component of 
the total UHECR flux becomes larger than the Galactic 
component should be somewhere around the so-called 
ankle, i.e. ~ 5 • 10 18 eV, while this transition would be 
at a lower energy in the pure proton case 0[3]|4]|5l. 

Any accurate description of the UHECR data should 
reproduce both the measured spectrum and composi- 
tion in a consistent way, and, while the measurement 
of the UHECR composition at Earth remains a difficult 
observational task, important progress has been made 
in the recent years. Notable results include those of 
the Pierre Auger Observatory |0 |7), which has pro- 
vided hints that the composition becomes heavier and 
heavier above ~ 10 19 eV (assuming the general valid- 
ity of hadronic interaction models). In contrast to this, 



other experiments such as HiRes (El and Telescope Ar- 
ray J9] ITOi have shown results which are compatible 
with pure proton scenarios. 

Propagation effects are known to modify the compo- 
sition of UHECRs, as the energetic nuclei are photo- 
dissociated in interactions with background photons. 
Horizon analyses show that Fe nuclei and protons with 
energies above ~ 6 ■ 10 19 eV can propagate over roughly 
the same distance without losing a significant fraction 
of their total energy, while intermediate mass nuclei are 
suppressed at shorter distances IfTTl fn. As a result UHE- 
CRs, at these high energies, should be dominated by ei- 
ther protons, Fe (or sub-Fe) nuclei, or a combination 
of the two. Therefore, the Auger results on UHECR 
composition can be understood if the proton component 
is cut at the source at a relatively low energy, around 
10 19 eV, and Fe nuclei are accelerated up to higher ener- 
gies, eventually dominating the overall spectrum. This 
would be natural, for instance, in a scenario where the 
different nuclei reach a maximum energy at the source 
that is proportional to their charge, Z. 

Although it has been shown that it is indeed possible 
to fit the UHECR energy spectrum in such a scenario 
0], this requires a source composition which is richer 
in Fe nuclei than would be expected from a simple ex- 
trapolation of the low energy cosmic ray source com- 
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position. This can be achieved, in principle, by invok- 
ing a source environment which is richer in Fe nuclei, 
or an acceleration mechanism which somehow discrim- 
inates nuclei on the basis of mass. In this paper, we pro- 
pose another mechanism to produce a heavier effective 
composition at the source, independent of the accelera- 
tion model. We assume only that the maximum energy 
reached by the particles accelerated in a given source is 
not universal and that the distribution of sources with 
respect to their maximum energy follows a power-law. 
The contribution of all sources is then equivalent to a 
scenario where identical sources not only inject UHE- 
CRs with a softer spectrum, as already shown by fl2l . 
but also with a heavier source composition. We then 
compute the Fe-to-proton enhancement ratio, rj Fe , as a 
function of the source parameters. 

2. Source E max Distributions and Resulting Spectra 

For most proposed cosmic ray sources the maximum 
energy which can be reached is limited by the ability of 
the source to contain particles in the acceleration region, 
i.e. by particle rigidity, defined as R = p/q. In the rela- 
tivistic limit applicable to UHECRs, this is proportional 
to E/Z, where Z is the charge of the accelerated nu- 
clei. In such a case, unless other mechanisms come into 
play to limit the energy of specific nuclei (such as photo- 
dissociation processes), the maximum energy, E^^, of 
nuclei of type i at the source is simply proportional to 
their charge, Z, : 
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where £^ is the proton maximum energy. 

For simplicity, UHECR models usually assume that 
all cosmic ray sources are identical, having the same 
spectrum extending up to the same energy. However, 
it is clear that the maximum energy will differ among 
sources depending on their individual properties, no- 
tably their size, magnetic field strength, intrinsic power, 
age, etc. In order to avoid introducing extra free param- 
eters the simplest assumption is to assume a power-law 
distribution for the number of sources as a function of 
E%i (as in 021): 
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where we have introduced E sup as the highest possible 
proton energy in any source and H{x) is the Heaviside 
step function. The parameter Eq is an arbitrary refer- 
ence energy, which has been introduced here simply to 



clarify the dimensionality of the various quantities. p 
is expected to be a positive number, as the number of 
sources able to accelerate nuclei to high energies pre- 
sumably decreases with increasing energy. 

We assume that the individual sources each produce a 
power-law spectrum of UHECRs with the same spectral 
index, x. The number of nuclei of type i injected per 
second and per unit energy by a given source is then: 
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where a,- is the abundance of nuclear species i, with 
I>; = 1. 

Assuming that the sources are homogeneously dis- 
tributed over space, the total number of cosmic rays in- 
jected per unit time, per unit energy, and per unit volume 
is then given by: 
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Replacing from Eqs.Q and Q, one gets: 



qi(E) = Q n a 
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which integrates (for f3 + 1) into: 
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NB: if p — 1, there are the same number of sources 
in each decade of E max , leading to a logarithmic diver- 
gence in the integral if E sup tends to infinity. In most 
natural cases, however, one would expect p to be larger 
than 1. 



3. Effective Spectrum and Composition in the limit 



of E « E. 
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For values of E suv much larger than E (and p > 1), 
Eq. (|6]l is equivalent to a single source power-law distri- 
bution given by: 



qi(E) = q ctjZf 
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In other words, the simple model described here is 
equivalent to the usual "universal source model" with 
the effective source spectral index 



x e ff = x + p — 1 , 



(8) 
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and with a modified source composition corresponding 
to the effective nuclear abundances: 



«i,eff = Ui X 4 
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The effect of a distribution of E lmx values among 
sources is thus a modification of both the energy spectral 
index and the abundances of nuclei in a correlated way 
through the parameter ft. This shows that these two as- 
pects of UHECR phenomenology are not independent. 

This behavior is easily understood by considering the 
simple example, say, of 5 discrete sources, as illustrated 
in Fig. [T] In this plot, the plain lines show an identical 
proton injection spectrum, extending up to a maximum 
energy that is different for each source. The dashed 
lines show the corresponding helium injection spectra, 
with an assumed abundance ratio of aH e /a p = 0.5, and 
a maximum energy twice as large for each source, as 
per Eq. ([T]). The bold lines are then simply the sum 
of all 5 contributions, showing an enhancement of the 
relative helium abundance. In this illustrative example, 
while the He nuclei are less abundant than the protons 
in each individual source, they dominate the overall flux 
injected at high energy. 

Fig. [T] also shows how the original source spectrum 
steepens as fewer and fewer sources contribute to the 
flux at higher and higher energies. Below the E mRX of the 
least energetic source, the sum spectrum obviously has 
the same form as the individual source spectra. Above 
that energy, and up to ~ E sup (i.e. the £ mal of the 
most energetic source), the resulting spectrum is a new 
power-law with the larger spectral index, x & g, as shown 
by the dashed line fit. 

In an actual astrophysical context, extensive propa- 
gation studies have determined the effective spectral in- 
dex needed to reproduce the observed UHECR energy 
spectrum for a given assumed source composition (see 
Sect.[T]i. In the case of pure proton sources the compo- 
sition effect demonstrated above is obviously irrelevant, 
but for a mixed composition scenario Eq. |9]l shows 
that the effective composition to be used in the mod- 
els is systematically richer in heavy nuclei than the in- 
dividual source composition, as soon as ft > 1. Given 
Eq. (|8j, this condition is equivalent to saying that the ef- 
fective source spectrum is steeper than the intrinsic one 
(x e ff > x). For instance, if the source spectral index is 
x = 2.0, as expected from standard diffusive shock ac- 
celeration, the effective source spectrum needed to re- 
produce the data in the case of a mixed composition 
model, namely x e g 2.3, requires a source distribution 
index ft - 1.3. This in turn implies that the effective 
Fe nuclei abundance is larger than the Fe abundance in 
individual sources by a factor Z^ 1 ~ 2.7. 
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Figure 1 : An illustration of the effect of an £ max distribution on the ra- 
tio of different nuclei in UHECRs. Here, we consider an artificial dis- 
tribution of sources with 5 different maximum energies, as explained 
in the text. The solid lines are the proton spectra for each source, while 
the dashed lines denote the He spectra, extending to twice higher en- 
ergies. The total effective spectrum for each species (before propa- 
gation) is given by the bold lines. In this example, even though the 
abundance of protons at the source is higher, in the UHECR region 
the spectra of the heavier nuclei are enhanced due to the distribution 
of the number of sources with respect to £ max . 



It is worth stressing that, as indicated above, the spec- 
trum given by Eq. (|7]i is formally that of a single source 
power-law distribution. Therefore, the propagated spec- 
trum that would be obtained from such a model, af- 
ter taking into account the interaction of the UHECRs 
with the intergalactic background fields, would fit the 
observed spectrum equally well, provided that the effec- 
tive spectral index (Eq. [8]) is adjusted in the same way 
as in the standard (single E max ) model. 

4. The Abundance Enhancement Ratio at the High- 
est Energies 

At the very highest energies, when E becomes closer 
to E sup , the last factor in Eq. (|6jl is no longer equal to 
~ 1 . The effective spectrum is then no longer a power- 
law, and the effective composition becomes dependent 
on energy. We can define the enhancement ratio of the 
abundance of the various nuclei relative to protons, rj,, 
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where q-,{E) was taken from Eq. ([6]). 

As an illustration, the behavior of the Fe enhance- 
ment ratio, 77Fe, as a function of E/E sup is plotted in 
Fig. [2] for several values of ft. The low energy limit is 
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Figure 2: The enhancement ratio of Fe nuclei, rjp e , Eq. jlO) , as a 
function of E/E sup , for several values of /3, the slope of the £ milx dis- 
tribution. 
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Figure 3: The enhancement ratio, 17;, for various nuclei, as indicated, 
as a function of E/E m „. 



given by 77, = zf 1 (see Eq. (|9jl). The enhancement ra- 
tio increases with energy. Likewise, the enhancement 
ratio is larger for heavier nuclei, as shown in Fig. [5] for 
,3=1.5. 

At the highest energies, close to iisup, the values 
shown on the plot should no longer be taken seriously, 
as in practice the granularity of the sources would have 
to be taken into account. The features of the cosmic ray 
energy spectrum and composition in this energy range 
are determined by the properties of the few highest en- 
ergy sources, which are subject to "cosmic variance". It 
is obvious, however, that at energies above the proton 
maximum energy, E sup , the UHECRs should be com- 
pletely devoid of protons, and thus the propagated UHE- 
CRs will be dominated by Fe and sub-Fe nuclei (TJ. 



5. Conclusion 

We have studied the effect of a distribution of 
UHECR sources as a function of the maximum en- 
ergy to which they can accelerate particles, both on 
the effective source spectrum and the effective source 
composition. Although the results shown here corre- 
spond to a specific set of assumptions, they show that 
the UHECR composition can indeed be significantly af- 
fected by such a distribution. These effects should thus 
be kept in mind when studying the phenomenology of 
UHECRs within the simplified framework of identical 
source models. 

In particular, we showed that a heavier composition 
can be obtained very naturally under the astrophysically 
sensible assumption that the number of sources able to 
accelerate particles up to a given energy decreases suf- 
ficiently rapidly with that energy (J3 > 1). The abun- 
dance enhancement effect depends on the charge of the 
nuclei, and is less pronounced for lighter nuclei. This 
tends to enhance the contribution of Fe nuclei with re- 
spect to intermediate mass nuclei, which are also more 
affected by propagation effects due to their shorter hori- 
zon scales m im . 

Such a heavier composition is needed to fit the most 
recent UHECR data, within a model where a low value 
of E sup is invoked to explain the apparent transition 
towards a heavy-nuclei-dominated component above a 
few 10 19 eV. These models naturally account for a rel- 
atively sharp cut-off of the proton-dominated UHECRs, 
in the energy range where the sources cease to acceler- 
ate protons. However, if the heavier components are not 
sufficiently abundant in this energy range, this proton 
cut-off would show in the energy spectrum as a visible 
feature. Therefore, the consistency of such models re- 
quires a higher abundance of heavy nuclei than the typ- 
ical expectations for the actual source composition. 

As we have shown here, an enhancement ratio of the 
order of 3 for Fe nuclei appears natural if the E max distri- 
bution function has a power law index ~ 1.3. This value 
of p is needed to go from an intrinsic source spectrum of 
E~ to an effective source spectrum of E~ 23 , as is typ- 
ically needed to fit the observed spectrum in the case of 
a mixed composition scenario 131 1131 . Larger enhance- 
ment ratios are possible if the E m . dx distribution func- 
tion is steeper, which then implies that the individual 
source spectra are harder. A source spectrum of 2? -1 - 8 , 
for instance, would need a value of /3 = 1.5 to mimic a 
single-type source distribution with a spectrum of E~ 23 , 
resulting in an enhancement ratio of ~ 5 for Fe nuclei. 

It should also be noted that, in the case of a low 
value of E sup , UHECR experiments are detecting cos- 
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mic rays just below and above E sap . The composition ef- 
fects would then be expected to be larger than discussed 
above (as apparent in Figs. [2] and |3J. In this transition 
range, across E snp , the analytical treatment used here is 
no longer relevant, and the local distribution of sources 
should be taken into account. For the same reason, a 
simple description of the overall spectrum in terms of 
power laws may not be possible for the highest energy 
cosmic rays. Even in the framework of the above con- 
tinuous model, Eq. (|6| is no longer a power-law in this 
energy range, and one cannot define an effective power- 
law index as in Eq. ([8}. 

We did not attempt a complete fit of the data here, as 
it will require additional assumptions about the individ- 
ual source composition(s) and possibly some fine tun- 
ing of the parameters of the most nearby sources. How- 
ever, in future works, we shall explore the cosmic vari- 
ance around E sup and investigate how the simple ideas 
studied here analytically can be applied to direct Monte- 
Carlo simulations with random source distributions, in 
order to create a global fit of present UHECR data, in- 
cluding both spectrum and composition. 
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